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Abstract. We study the impact of 1// noise on the 
Planck Low Frequency Instrument (LFI) observations 
(Mandolesi et al. 1998) and describe a simple method for 
removing striping effects from the maps for a number 
of different scanning strategies. A configuration with an 
angle between telescope optical axis and spin-axis just 
less than 90° (namely ~ 85°) shows good destriping ef- 
ficiency for all receivers in the focal plane, with resid- 
ual noise degradation < 1 — 2%. In this configuration, 
the full sky coverage can be achieved for each channel 
separately with a 5° spin-axis precession to maintain a 
constant solar aspect angle. 

Key words: Methods: data analysis - Cosmology: cos- 
mic microwave background - space vehicles 
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1. Introduction 

The great success of COBE-DMR results (Smoot 
et al. 1992, Bennet et al. 1996a, Gorski et al. 1996) had 
a considerable impact on the scientific community being 
the first detection of CMB (Cosmic Microwave Back- 
ground) anisotropy and supporting the gravitational in- 
stability scenario for structure formation. Detailed pre- 
dictions for the CMB angular power spectrum in dif- 
ferent cosmological models and the proper signature of 
different cosmological parameters (e.g. Zaldarriaga & 
Seljak 1997) are now quite well established. A signif- 
icant number of new experiments, both ground-based 
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and balloon-borne, with sensitivity and angular reso- 
lution better than COBE (e.g. Lasenby et al. 1998 and 
De Bernardis & Masi 1998) are going to be designed and 
constructed and will deliver good quality data sets in the 
next few years. 

It is also clear (Danese et al. 1996) that only a space 
mission free from unwanted contamination from ground 
and Earth atmosphere and with a nearly full-sky cov- 
erage can fully exploit the gold mine of cosmological 
information imprinted in the CMB anisotropy. The two 
space missions MAP (Microwave Anisotropy Probe) (see 
Bennet et al. 1996b) by NASA and Planck (Mandolesi 
et al. 1998, Puget et al. 1998) by ESA are planned for 
the next decade, the first to be launched in the year 
2001 and the second in 2007. These planned observa- 
tions, with their wealth of cosmological information, will 
play a determinant role to fully understand the prop- 
erties and the evolution of the universe. In particular, 
Planck is a third-generation CMB space missions de- 
signed to achieve unprecedented sensitivity and angular 
resolution over the entire sky. The extremely wide fre- 
quency range (~ 30 -^ 900 GHz) explored by the com- 
bination of two focal plane arrays, the Low Frequency 
Instrument (LFI) and the High Frequency Instrument 
(HFI), will allow to efficiently separate the frequency in- 
dependent primary cosmological information from astro- 
physical signatures (e.g. De Zotti et al. 1999) by taking 
advantage from their proper frequency dependencies. 

As for any CMB experiment, great attention has to 
be devoted by the two mission teams to all the possible 
systematic effects. In the context of the Planck mis- 
sion we have carried out a detailed study of one of these 
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effects: the so-called 1//, or low-frequency, noise which 
is particularly relevant for the LFI radiometers, lead- 
ing to stripes in the final maps which increase the over- 
all noise level and may alter the statistical analysis of 
the anisotropy distribution. Recent analytical work by 
Seiffert et al. (1997) has shown the dependency of the 
1// noise upon the radiometer characteristics such as 
the bandwidth, the noise temperature, payload environ- 
ment temperature and other quantities properly related 
to Planck-LFI radiometers. They can be combined to 
define a representative parameter, the "knee-frequency" 
fk, which has to be kept as low as possible compared 
with the spinning frequency fs of the spacecraft. Janssen 
et al. (1996) have indeed demonstrated that for fk ^ fs 
a degradation in final sensitivity will result. 

In Sect, y we briefly summarize the source of the 
1// noise and the dependency of the knee-frequency 
upon the radiometer properties. In Sect. ^ we outline 
the framework of our simulations, evaluate the possible 
impact of 1// noise and present the technique for re- 
moving the 1// noise in the context of Planck-LFI. 
The main results both for un-reduced and reduced 1// 
noise are presented in Sect. 0. Our main conclusions and 
implications for the scanning strategy are discussed in 
Sect. |. 

2. Source of i// noise 

The presence of unwanted systematic effects is one of the 
main difficulties in CMB anisotropy experiments. The 
so-called 1// noise, typically generated by HEMT (High 
Electron Mobility Transistor) amplifier gain instabilities, 
is of particular interest for the Planck-LFI since, if not 
properly corrected, it may lead to stripes in the final 
maps due to the satellite scanning strategy. Therefore it 
is of great importance to reduce the impact of such ef- 
fect both by hardware and software techniques. The LFI 
receiver concept is driven by the need to reduce instabil- 
ity effects. Bersanelli et al. (1995) described the design of 
LFI radiometers (Fig. yj) which are modified Blum corre- 
lation receivers (Blum 1959, Colvin 1961) known also as 
pseudo-correlation or continuous comparison receivers. 
The details of the radiometer will not be described here. 
The main feature to note is that the receivers have two 
inputs: one from one of the two polarization directions of 
a feed-horn and one from a reference load. The signal is 
sent by an hybrid coupler in the two "legs" of the receiver 
where it is amplified and phase shifted. The signal after a 
second hybrid, is amplified again and detected. The out- 
put signal is proportional to the signal difference at the 
two inputs. This is achieved by differencing successive 
measurements at the detector which, for the receiver de- 
sign, come from the horn and from the load respectively. 
The ideal situation with respect to the sensitivity to sys- 
tematic effects is when the output is nulled. It is there- 



Feed Horn 



Hybrid Coupler 



Phase Switch 



HEMT Amplifiers 




Reference Load 



InPHEMT 



Fig. 1. Schematic representation of the LFI radiometers. 

fore important that the temperatures of the reference 
load be as close as possible to the « 3 K temperature 
of the sky signal. Furthermore, a nearly-null output is 
obtained by adjusting properly the DC gain ratio r after 
the diodes at the value r — [T^ -t- Tn)/{Ty + Tn), where 
Ty is the temperature from the reference load input. 

Seiffert et al. (1997) have reported an accurate anal- 
ysis of the source of 1// noise in the final output of this 
receiver design. The ideal white noise sensitivity for this 
kind of radiometer is given by: 



AT — 



\/2(r„ + T,) 



/Bt 



(1) 



where B is the bandwidth (for LFI 20% of the central 
frequency), r is the integration time, T^ is the signal at 
horn input and r„ is the noise temperature of the ampli- 
fier. The main source of 1// noise resulted be gain fluc- 
tuations in the HEMT amplifiers which have a typical 
1// spectrum. These fiuctuations translate into fluctua- 
tions of the noise temperature T„ with the same spectral 
shape. Noise temperature fiuctuations Ar„ can mimic a 
change in the radiometer output 
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where A is the square root of the spectral ampli- 
tude of noise temperature fluctuations. Following Seif- 
fert et al. (1997), we define the "knee-frequency" when 
This leads to the expression 
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Typical values are ^ -- 1.8 x 10"'"' for 30 and 44 GHz 
radiometers and A -^ 2.5 x 10"^ for 70 and 100 GHz. 
Fig. g reports an example of simulated output for a pure 
white noise and a white noise plus 1// noise signal: their 
difference is much more clear in Fourier space where the 
1// noise spectrum shows up at low frequencies. For the 
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Fig. 2. Comparison between a pure white noise (right) and white+1// noise (left) in time domain (upper panels) and and 
Fourier space (bottom panels). Their difference is more clear in Fourier space with a flat spectrum for the white noise and 
the typical 1// noise spectrum for the low-frequency noise. We used a, fk = 0.8 Hz for this illustrative case 



30 GHz channel, with T„ = 9 K, ^^ = 3 K, 5=6 GHz, 
Ty ^ 20 K Eq.d) gives a /fe ~ 0.046 Hz. For a 4 K 
reference load, which is now the adopted basehne, the 
knee-frequency reduces to « 1 niHz. For the 100 GHz 
channel fk is higher (0.11 Hz) for a 20 K reference load 
and is reduced to less than 1 mHz for the 4 K load. 



3. Scan circles, sky maps and destriping 

As reported in Mandolesi et al. (1998), the selected or- 
bit for Planck satellite will be a Lissajous orbit around 
the L2 Lagrangian point of the Sun-Earth system. The 
spacecraft spins at 1 r.p.m. and the spin axis is kept 
on the Ecliptic plane at constant solar angle by repoint- 
ing of 2.5' every hour. The field of view of the two in- 
struments is between a w 80° — 90° from the spin-axis 
direction. Hence Planck will trace large circles in the 
sky: these circles cross each other in regions close to the 
Ecliptic poles. The shape and width of these regions de- 
pend upon the angle a, the scanning strategy and beam 
location in the focal plane. The value of the angle a has 
not yet been fully defined, as well as the scanning strat- 
egy, which may or may not include a periodic motion 
of the spin-axis away from the Ecliptic plane. These op- 
tions depend on a trade-off between different systematic 
effects (striping, thermal effects, stray light), which has 
to be carefully addressed. 

Here we made an attempt in this direction taking 
into consideration the 1// noise alone as if all the other 
systcmatics were negligible. This is of course an idealised 
situation but gives us an estimate of the impact of this 



effect that therefore should be regarded as a lower limit. 
Nonetheless this result may help for a better design of 
the instrument. Future work will consider a more com- 
plete situation when a detailed thermal transient analy- 
sis will be available for the Planck spacecraft. 

3.1. The "flight- simulator" 

Burigana et al. (1997, 1998) has described in detail the 
code we have implemented for the Planck scanning 
strategy and we refer the reader to these papers. The 
relevant geometrical inputs are the beam location in the 
focal plane and the angle a between the spin and point- 
ing axis. 

For each beam position on the focal plane our code 
outputs the complete data stream. We consider here a 
reduced version of the actual baseline for the scanning 
strategy (actual parameters in parentheses): spin-axis 
shift of 5' every 2 hours (instead of 2.5' every hour) and 
three samplings per FWHM of 30' at 30 GHz (instead 
of 12 samplings every 30', i.e. 4 samplings every 10', the 
FWHM at 100 GHz; see Mandolesi et al. 1998). These 
modifications allow us to explore a large region of the 
parameter space, beam position, fk, scanning strategy 
and pointing angle a, in reasonable time. Furthermore 
we do not consider the single minute data stream but 
we take the average over the 120 circles forming a given 
2-hours set. In what follows we run simulations for the 
30 GHz channel. 

Wright (1996) has shown that possible data filtering 
on a given scan circle may help in reducing the impact 
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of 1// noise. This is useful for values of fk typical for 
"total power" receivers which are much higher that our 
0.05 Hz and we chose not to include this technique here. 
It is clear from Eqs.(|l|) and (^) that both the sen- 
sitivity and the "knee-frequency" depend on the actual 
temperature in the sky T^ seen by the horn. Our syn- 
thetic model for microwave sky emission includes a stan- 
dard CDM prediction for CMB fluctuations plus a model 
of galactic emission. This model has the spatial tem- 
plate from the dust emission (Schlegel et al. 1998) but 
has been normalized to include contribution from syn- 
chrotron, free-free and dust according to COBE-DMR 
results (Kogut et al. 1996). The major foreground con- 
tamination at 30 GHz comes from synchrotron and free- 
free. We then choose to overestimate the overall syn- 
chrotron fluctuations by a factor of ~ 10. This is the 
worst case scenario with respect to destriping efficiency 
(see Sect. O and O). 



Of course the impact of the sky temperature is small, 
since T^ is a small fraction of the noise temperature r„ : 
this is in fact ~ 10 K while Tj, « 3 K at 30 GHz and the 
maximum galactic contamination adopted here is about 
44 mK. In this evaluation of T^ we include also a typical 
environment temperature of about 1 K. 

We do not perform a convolution with a real beam 
pattern and sky signal since the real beam pattern will be 
in general not perfectly symmetric, not perfectly gaus- 
sian and with sidelobe contributions which will introduce 
other kinds of contamination. Since we want to estimate 
here only the impact of 1// noise we convolve our input 
map with a gaussian beam with the nominal FWHM 
(33') of the 30 GHz Planck-LFI channel we consider. 

3.2. Generation of instrumental noise 

We have the possibility to generate different kinds of 
noise spectra. We work in Fourier space and generate 
the real and imaginary part of Fourier coefficients of our 
noise signal. The noise spectrum has the usual form: 
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where a is a normalization factor and /? takes typical val- 
ues from 1 to 2.5 depending on the source of noise (gain 
drift of thermal effects): /3 = 1 is the 1// noise case. In 
general it is a good approximation to take Snoiseif) van- 
ishing for / < fjni7i and / > f,nax- Delabrouille (1998) 

proposed Jrain ^^ ^/ -^mission anCL Jmax ^^ ^ / -^-^ samplings- 

After generating a realisation of the real and imag- 
inary part of the Fourier coefficients with spectrum de- 
fined in Eq.(|j), we FFT (Fast Fourier Transform; Cooley 
& Tukey 1965; Heideman et al. 1984) them and obtain a 
real noise stream which has to be normalized to the white 
noise level defined in Eq.dl^). We use a typical value of 
fk = 0.05 Hz for the knce-frcquency at 30 GHz assuming 



a 20 K load temperature (see Sect. 0). We chose to gen- 
erate noise streams of about 2 x 10 data points, corre- 
sponding to 8 spin axis positions or 16 hours: this seems 
a reasonable compromise among our present knowledge 
of real hardware behaviour and the required computa- 
tional time with respect noise stream length. The actual 
time for a generation of a noise stream 16 hours long 
and one year of mission even in our reduced scanning 
strategy will require about 10 hours on a Silicon Graph- 
ics machine with 2 Gb of RAM and clock speed of 225 
MHz. We also verified that the most time consuming op- 
eration in our code is just the FFT 1// noise generation. 
The possibility to generate the noise streams with a 
completely arbitrary spectrum in real space (e.g. Becca- 
ria et al. 1996 and Cuoco & Curci 1997 and references 
therein), without the use of FFT, seems very promising 
for considerably reducing the computational time (Wan- 
delt et al. 1999, Wandelt & Gorski 1999). 

3.3. From data stream to sky maps 

The final output of our "flight-simulator" are 4 matrices 
with a number of rows equal to the considered spin-axis 
positions ris (for one year of mission ~ 4320 in our re- 
duced baseline as in Sect.pj) and a number of columns 
equal to the number of integrations, weakly dependent 
on a, along one scan circle (here Up ~ 2160). The input 
and output maps are in HEALPix pixelisation sch eme 
(Gorski et al. 1998, http: //www.tac .dk/'healpix ); we 
use an input map with a resolution of about 3.5' cor- 
responding to 12 million equal area pixels on the full 
sky. 

For each circle, the code outputs are the pixel num- 
ber N at the specified resolutions, the temperature plus 
total noise contribution T (white plus 1// noise), the 
temperature with only white noise W and the pure sig- 
nal as observed in absence of instrumental noise G. The 
W and G will be used for studying the degradation of 
1// noise with respect to the ideal pure white noise case 
and the impact of scanning strategy geometry on ob- 
served pixel temperatures. 

We can arbitrarily choose the temperature output 
data stream resolution from 3.5' to higher values (smaller 
resolution) which set also the output temperature map 
resolution. Regarding the data stream for the pixel num- 
ber outputs we can also use higher resolutions, allowing 
to test the impact of using more or less stringent crossing 
conditions in the destriping algorithm (see the following 
section). 

From these data streams it is quite simple to ob- 
tain observed simulated maps: we make use of N and 
T to coadd the temperatures of those pixels observed 
several times during the mission. In the same way we 
build maps with only white noise contribution, without 
receiver noise, as well as a sensitivity map knowing how 
many times a single pixel is observed. 
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Fig. 3. Pure noise (white and 1// noise together) map before destriping: stripes are clearly present. The adopted fk is 0.05 
Hz and units are in mK. See the text for simulation input parameters. 

In Fig. ^we show a pure noise (white plus 1//) map the result will be matrices with Ug x ni rows and Up/ni 
in Ecliptic coordinates after signal subtraction (T — G): columns, 
stripes are clearly present. 



3.4- Destriping techniques 

We developed a simple technique which is able to elimi- 
nate gain drifts due to 1// noise. This is derived from the 
COBRAS/SAMBA Phase A study proposal (BersaneUi 
et al. 1996) and from a re-analysis of Delabrouille (1998). 
As reported by Delabrouille (1998) the effect of 1// noise 
can be seen as one or more additive levels, different for 
each scan circle. We worked with averaged (over 2 hours 
period) scan circles and hence we nearly removed drifts 
within each circle: what is left is related to the "mean" 
1// noise level for this observation period. In fact aver- 
aging scan circles into a single ring corresponds to a low- 
pass filtering operation. As long as fk is not far larger 
than the spin frequency, this ensures that only the very 
lowest frequency components of the 1// noise survive. 
Therefore it is a good approximation to model the aver- 
aged 1// noise as a constant offset Ai for each ring for 
the set of parameters we are using. 

We want to obtain the baselines for all the circles 
and re-adjust the signals correspondingly. Delabrouille 
(1998) claimed the possibility that working with multiple 
baseline per circle (n; — number of baseline levels) may 
result in a better destriping. In our code this is equivalent 
to dividing the rows of each output matrices in ni parts: 



In order to estimate the different Ai we have to find 
common pixels observed by different scan circles and the 
pixel size in the matrix N is a key parameter. Increasing 
the resolution used in N reduces the number of cross- 
ings possibly yielding to lower destriping efficiency, while 
adopting resolutions lower than the resolution of the in- 
put map and of the matrix T introduces extra noise, 
related to variations of real sky temperature within the 
scale corresponding to the lower resolution adopted in 
N, which is not taken into account in the estimate of re- 
ceiver noise and may introduce artifacts in the destriping 
code. The adopted pessimistic galactic emission model, 
which by construction has gradients larger than those 
inferred by current data, emphasizes this effect and our 
simulations are then conservative in this respect. In the 
following Nil, Til and En will denote the pixel number, 
the temperature and the white noise level for the pixel 
in the i**^ row and l^^ column. Let us denote a generic 
pair of different observations of the same pixel with an 
index tt which will range between 1 and Uc, the total 
number of pairs found. In this notation it is related to 
two elements of N: tt — > {il,jm) where i and j identify 
different scan circles and I and m the respective position 
in each of the two circles. 
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We want to minimize the quantity: 
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with respect to the unknown levels Ai. The sub- index 
TT indicates that each set of {il,jm) is used in that 
summation. It is clear that S is quadratic in all the 
unknown Ai and that only differences between Ai en- 
ter into Eq.(5). Therefore the solution is determined 
up to an arbitrary additive constant (with no physi- 
cal meaning for anisotropy measurements). We choose 
then to remove this indetermination by requiring that 
J2h=i ^h = 0- This is equivalent to replace Eq.(5) with 

S' - 
get: 



S + V^hJlAh) ■ After some algebra we finally 
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for all the k = 1, ...,ns (here 5 is the usual Kronecker 
symbol) . This translates into a set of n^ linear equation 
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which can be easily solved. Furthermore we note that by 
construction the matrix C of Ckh coefficients is symmet- 
ric, positive defined and non singular. The first property 
permits to hold in memory only one half of the matrix 
C (e.g. the upper- right part), the second allows us to 
solve the linear system without having to exchange rows 
or columns (Strang 1976), so preserving the symmetry. 
The non-singularity of C is true provided that there are 
enough intersections between different circles and hence 
is related to the resolution at which we look for com- 
mon pixels between different scan circles, to the scan- 
ning strategy and beam location. A detailed discussion 
of numerical algorithm for solving this system with sig- 
nificant saving of required RAM is presented in Burigana 
etal. (1997). 

It is interesting to note that the applicability of this 
destriping technique does not depend upon any a-priori 
assumption about the real value of /^ or the real noise 
spectral shape since it can work also for different val- 
ues of the exponent /? in Eq.(^), althought its efficiency 
could be in principle partially optimized in relation to 
the properties of the noise and of the sky gradients. 




Fig. 4. Noise power spectra at 30GHz before destriping. Sim- 
ulations parameters are: a = 90°, {9b, (^b) = (2.8°, 45°), 
fk = 0.05 Hz and spin-axis always on the Ecliptic plane. 
Theoretical white noise level is also reported for comparison. 
The excess of noise is about 40% over the white noise level. 



(6) 4. Simulations results and scanning strategy 



We first consider an angle a = 90° between telescope 
and spin axes and a beam location with (9b,4>b) = 
(2.8°, 45°). Here 9b is the angle from the optical axis: 
this is a typical value for the 100 GHz (with a ~ 1.5 m 
aperture off-axis Gregorian telescope) horns while the 
30 GHz beams are placed at larger Ob values. This as- 
sumption is therefore conservative with respect to the 
destriping efficiency, since in this case the region of cross- 
ings between scan circles is smaller and closer to Eclip- 
tic poles. The angle (J)b is the beam center displacement 
from the axis given by the intersection between the sky 
field of view plane and the plane containing the telescope 
and spin axes (see Burigana et al. 1998). Our choice of 
(l)B = 45° is intermediate between 0° (or 180°) and 90° 
which are equivalent for the destriping respectively to 
an on-axis beam (with crossings only at Ecliptic poles 
for a = 90°) and a beam which spreads crossings over 
the wider possible region. The adopted beam location is 
therefore a non-degenerate case although non optimal. 
The knee frequency fk = 0.05 Hz and the microwave 



sky emission is modeled as in Sect. 3.1 



We evaluated the impact of 1// noise both in terms 
of added rms noise and of angular power spectrum. 
The output of our simulations is expressed in antenna 
temperature and has to be converted into thermody- 
namic temperature for the comparison with CMB power 
spectrum, as usually predicted in theoretical models. 
At 30 GHz the conversion factor is very close to unity 
(rmsth — rms A, where rms a is the rms in antenna 
temperature). 
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Fig. 5. Baselines recovered from destriping algorithm for the case reported in Fig 



The white noise power spectrum can be derived an- 
alytically knowing the total number of pixels in the sky 
and pixel sensitivity : 
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(8) 



J^.l. Morphology of the C^ artefacts due to 1// noise 

We now give a short theoretical argument for the ex- 
pected morphology of the Ci power spectra of the simu- 
lated striping pattern. For scanning without wobbling or 
precession, two angular scales fully determine the struc- 
ture of the striping pattern on the sky. One is the half 
opening angle Ai = 90° of the ecliptic equator which 
the spin axis traces on the sky as Planck completes a 
full orbit about the sun. This makes the noise pattern 
symmetric under parity, on average. It follows that the 
correlation function for the noise pattern is symmetric 
about 90°, with a strong peak at 180° and for the 1// 
component CV = for all odd L 

The other scale is the half opening angle of the scan- 
ning rings A2 — a — 6'_B sin(0B), 9b ^ 7r/2. This leads 
to peaks in the correlation function at 2A2 and by sym- 
metry at 180° - 2A2. 

The combination of these two nearby scales leads 
to beats in the C^, which are visually apparent as the 
"blob" in Figs. 4,8,10,12 and 14. Treating this effect in 
the same way as the appearance of fringes in an interfer- 
ence pattern we can calculate the width in t of each blob 



as A^ = 90/ (A2 — Al). This is in quantitative agreement 
with our numerically computed spectra. It also explains 
the absence of beats in Figure 10, where A2 = Ai = 90° 
and the power spectrum can be understood as a single 
blob with infinite /S.i. 

For scanning strategies with some kind of "wobble" , 
ie. where the spin axis does not trace out the ecliptic 
equator on the sky, the correlation function will still be 
approximately symmetric if the amplitude of the wob- 
ble is small compared to 90° and the above statements 
remain approximately true. However, the fact that the 
symmetry is weakly broken can be seen in Figure 8, 
where the interference effect is slightly washed out and 
the beats do not reach all the way down to the white 
noise level. 

The offset of the blobs compared to the white noise 
level is given by the rms excess power due to the 1// 
noise component. At first sight it may be surprising 
that other than this normalisation, the shape of the 
excess power power spectrum Ci of the striping pat- 
tern is independent of the spectral content of the noise, 
as parametrised by the knee frequency. This is due to 
the filtering effect of co-adding scan circles into rings 
discussed above. In the regime which our simulations 
probe, the structure which remains on each ring is well 
described by a constant offset plus white noise. Hence, 
as long as sufficiently many scan circles are co-added, 
changing the knee frequency simply scales the relative 
contribution of 1// and white noise to the noise pattern 
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Fig. 6. Noise map after destriping for the same case illustrated in Fig. ra stripes are no more visible. 



4-2. Simulation results 

In Fig. H we show the square root of the power spectrum 
(which is roughly proportional to rms contribution of 
temperature per ^-bin) before applying the destriping 
technique: the solid line is the white noise level as de- 
rived from Eq. (0) and the superimposed spectrum is de- 
rived from a simulation with only white noise included. 
The agreement is very good which confirms the accuracy 
of our simulator and map-making algorithm. The gray 
line is the global noise spectrum (white and 1// noise 
together). The spectrum is clearly non- white: blobs are 
present. The excess of noise in terms of both rmsth and 
y/U^ is about 40% of white noise level. 

Fig. |g shows the noise map after applying our de- 
striping code: stripes are no more evident and Fig. || 
shows the recovered baselines by destriping algorithm. 
This situation can be quantified by computing the noise 
power spectrum of the destripped map, as shown in 
Fig. 1^. Now the spectrum is considerably flatter and no 
more blobs are present, with an overall noise excess of 
1-2% over the white noise level. This confirms the ef- 
ficiency of the destriping algorithm, under the above 
mentioned simplifying assumptions. On the other hand 
it is clear from Fig. M that at large scales, correspond- 
ing to £ ;S 100, a non-negligible residual contribution is 
present in the noise spectrum. A lower knee frequency 
(achievable with the reference loads at 4K) is necessary 
to further reduce the effect. We also investigated if differ- 
ent scanning strategies may help in destriping efficiency. 



Possible periodic motion of the spin-axis away from the 
Ecliptic plane have the effect of broadening the region 
of crossings between different circles. This goes in the 
direction of removing possible degeneracies in the de- 
striping system. We implemented both sinusoidal oscil- 
lations and precessions: the first does not preserve the 
spacecraft solar illumination and this is likely to induce 
thermal effects and drifts in the data. 

The second is then preferable, since it keeps the solar 
angle constant as the satellite moves the spin-axis. For 
both the motions we performed 10 complete oscillations 
per year of mission with 10° amplitude. The results be- 
fore and after destriping are nearly the same for the two 
cases and Fig. @ and Fig. M report noise power spectra 
for the precession motion: before destriping blobs are 
still present and disappear after destriping. The excess 
of noise is ~ 40% and 1-2% before and after destriping 
respectively. 

We consider also a case of an on-axis beam 
{ObAb) = (0°, 0°) with a = 90°. This situation is repre- 
sentative of horns which, in the focal plane arrangement, 
are placed close to the scanning direction: from the de- 
striping point of view these horns are indeed equivalent 
to on-axis horns and intersections between different scan 
circles are nearly only at Ecliptic poles. This represents 
a degenerate situation with respect to destriping effi- 
ciency. Fig. no and Fig. 11 show noise power spectra for 
this on-axis case before and after destriping respectively: 
the geometry of the simulation is changed and the same 
for blob shape which is not well defined. Here the level of 
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Fig. 7. Noise power spectra at 30GHz after destriping. Sim- Fig. 8. Noise power spectra at 30GHz Before destriping. Sim- 
ulation parameters are the same of Fig. H. Now the added ulation parameters are the same of Fig. MJ except for the scan- 
noise is only 1-2% of the white noise level. White noise spec- ning strategy which includes here precession (see the text), 
tra are also reported for comparison. The added noise is ~ 40% of the white noise level. White 

noise spectra are also reported for comparison. 



added noise is ^ 42%. After destriping we are left with 
an excess of noise of the order of 5% of the white noise 
and the noise spectrum is considerable less flat than be- 
fore. Such excess noise and a non-flat spectrum around 
^ - 100 - 200 (where the first CMB Doppler peak is ex- 
pected) after destriping is clearly not acceptable for the 
feeds aligned with the scan direction. Moving the spin 
axis is a way to remove this degeneracy, as well as to 
complete the sky coverage for all channels. Keeping fixed 
the beam position we move the angle a from 90° to 85°: 
this of course leaves small regions around Ecliptic poles 
which are not observed but has the advantage to enlarge 
the scan circle crossing region. Fig. ^ and Fig. |l^ show 
noise power spectra before and after destriping respec- 
tively for this 85°, on-axis case: now blobs are clearly 
visible with dimensions in ^-space different from Fig. ^ 
due to different geometrical configuration. These blobs 
are completely removed after destriping leaving a level 
of added noise after destriping of only ;$ 2% over the 

white noise. This is essentially what we obtained for the 
off-axis configuration. All the simulations considered so 
far were performed with a fk — 0.05 Hz which is derived 
from Eq.(||) assuming a Tj, = 20 K of reference load. 
To decrease the residual effect, a 4 K reference load is 
being designed, exploiting the 4 K stage of the HFI. In 
this configuration it is expected that the resulting knee 
frequency would be less than 10 mHz. We have run a 
simulation with f^ ~ 0.01 Hz with the usual off-axis 
configuration with a = 90°. As for the other cases we 
report noise power spectra before and after destriping in 
Fig. nj and Fig. n5l It is interesting to note that now the 



level of added noise is 9% before destriping and reduces 
only to 0.5% after applying destriping algorithm. 

Although a comparison among different kinds of sys- 
tematic effects (Burigana et al. 1999) and studying how 
they combine in LFI observations is out of the aim of the 
present work, we have tried a first attempt in considering 
the sensitivity of this destriping algorithm to one of the 
most relevant source of contamination in Planck mea- 
surements, the beam shape asymmetry. We have con- 
sidered two simulations: one with a pure circular gaus- 
sian beam and another with a highly elliptical gaussian 
beam with a ratio of 1.7 between major and minor axis, 
that exibit rms differences of about 5 /zK in the tem- 
peratures observed in the same sky directions (Burigana 
et al. 1998). The impact of beam ellipticity on destrip- 
ing efficiency appears to be small: it has a maximum 
effect which is ~3 times smaller than the white noise 
level between 100 < £ < 200 (Burigana et al. 1999). This 
preliminary result suggests that the adopted algorithm 
may be only weakly dependent on the contamination 
introduced by other kinds of systematic effects; further 
investigations are demanded to confirm this perspective. 

5. Discussion and conclusion 

We have reported here an extensive study of the 1// 
noise contamination for different scanning strategies, 
beam location on the focal plane as well as for differ- 
ent values of the knee frequency. Even in this idealised 
situation (all other systematics are well under control) 
we can gain useful information from our study. 
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Fig. 9. Noise power spectra at 30GHz after destriping. Sim- 
ulation parameters are the same of Fig. H. The added noise is 
now only 1-2% of the white noise level. White noise spectra 
are also reported for comparison. 



First of all from these simulations it seems that mov- 
ing the spin-axis away from the ecliptic plane does not 
significantly help the destriping efficiency for typical LFI 
beam locations and, concerning the 1// noise alone as 
source of drifts, it would be preferable to keep the spin- 
axis always on the Ecliptic plane. This is clear when we 
compare results as those in Fig ^ and Fig. which are 
practically identical. 

Furthermore for most of the LFI beams the choice 
of a = 90° would be acceptable: the destriping left only 
;S 2% of excess noise with respect to the pure white noise 
case. On the other hand, some LFI beams are equivalent 
to on-axis beam which, from Fig. |l^, is clearly a degen- 
erate case. A smaller value of a (namely 85°) breaks this 
degeneracy at an acceptable level yielding the usual de- 
striping efficiency. From these two points an indication of 
a possible choice of the scanning strategy and instrument 
configuration arises: with a = 85° and precession of the 
spin- axis (no thermal drifts) with only 5° amplitude (half 
of what we considered here) appears satisfactory for de- 
striping performances while preserving full-sky coverage 
for all channels. This allows data redundancy but intro- 
duces irregularities in the integration time distribution, 
which may be an issue for the data analysis. Without 
spin axis modulations, a quite complete sky coverage 
and a smooth integration time distribution at each fre- 
quency can be achieved only by assembling data from 
different receivers, losing redundancy. 

We have considered also the destriping efficiency 
when using more than a single baseline level Ai per scan 
circle: results reported in Burigana et al. (1997) indicate 
that for the 1// noise case no significant improvements 



Fig. 10. Noise power spectra at 30GHz Before destriping. 
Simulation parameters are the same of Fig. Ml but the beam 
position is now {6b, (^b) = (0°, 0°). The added noise is ~ 42% 
of the white noise level. White noise spectra are also reported 
for comparison. 



are achieved: either the 1// noise level per circle is al- 
most constant and equivalent to a single baseline, or the 
accuracy in the determination of more levels per circle 
decreases. The whole set of simulations seems to indi- 
cate that there is enough redundancy of observations to 
remove at acceptable level the contamination due to 1// 
noise even if we require a more strict condition of cross- 
ings between scan circles. Of course, the performance of 
this destriping code could be partially optimized in the 
future by appropriately choosing the number of levels 
per circle and the crossing condition according to the 
dominant kind of instrumental noise (the parameters fk 
and /3), the magnitude of the gradients in the sky emis- 
sion and our knowledge of other contamination sources. 
For what concerns properly the 1// noise, an im- 
portant indication comes from the simulation with fk = 
0.01 Hz: the excess of noise before destriping reduces by 
a factor ~ 4 with respect to the case fk = 0.05 Hz. 
Reaching fk = 0.01 Hz is a goal of the LFI design, which 
will ensure a clean and robust observation. In addition 
the extra noise level after destriping decreases, at least 
under these ideal assumptions, by a factor ~ 3. 

There are many open issues both astrophysical and 
instrumental. Regarding the first, the microwave emis- 
sion model we use, although pessimistic for what con- 
cerns galactic emission, can be completed with the in- 
clusion of different foreground contributions. The emis- 
sion from extragalactic point sources and in particular 
their variability may decrease our destriping efficiency. 
We have also not considered here any other source of 
possible systematic effects such as thermal effects and 
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Fig. 11. Noise power spectra at 30GHz after destriping. Sim- 
ulation parameters are the same of Fig. hcl. The added noise 
is now larger than previous cases, being 5% of the white noise 
level. White noise spectra are also reported for comparison. 



stray-light contamination induced by Galaxy emission. 
Burigana et al. (1999) describe in detail the impact of 
stray-light contamination. Also these effects may in prin- 
ciple degrade the accuracy in removing 1// noise stripes. 
As proved by the scientific experience in many years 
of work both in physics and in cosmology and astro- 
physics, efficient data analysis tools can significantly im- 
prove the quality of the information extracted from the 
data, provided that the systematic effects are well un- 
derstood, but the first and most important step in pro- 
jecting experiments is to reduce all the contaminations 
at the lowest possible levels. It is then of great impor- 
tance to decrease as much as possible the impact of 1// 
noise before destriping and fk = 0.01 Hz is an important 
goal for instrument studies and prototypes. 
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Fig. 12. Noise power spectra at 30GHz Before destriping. 
Simulation parameters are the same of Fig. hd but the an- 
gle Q is now set to 85°. The added noise is ~ 42% of the 
white noise level. White noise spectra are also reported for 
comparison. 
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Fig. 13. Noise power spectra at 30GHz after destriping. Sim- Fig. 15. Noise power spectra at 30GHz after destriping. Sim- 
ulation parameters are the same of Fig. [13. The added noise ulation parameters are the same of Fig. fWl The added noise 
is now only 2% of the white noise level. White noise spectra is now only 0.5% larger the white noise level. White noise 
are also reported for comparison. spectra are also reported for comparison. 
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Fig. 14. Noise power spectra at 30GHz before destrip- 
ing. Simulation parameters are the same of Fig. H but now 
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